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Excess deaths provide total impact estimates of major crises, such as the COVID-19 
pandemic. We evaluated excess death trajectories across countries with accurate 
death registration and population age structure data and assessed relationships with 
vulnerability indicators. Using the Human Mortality Database on 34 countries, 
excess deaths were calculated for 2020—2023 (to week 29, 2023) using 2017-2019 
as reference, with adjustment for 5 age strata. Countries were divided into less and 
more vulnerable; the latter had per capita nominal GDP < $30,000, Gini > 0.35 
for income inequality and/or at least >2.5% of their population living in poverty. 
Excess deaths (as proportion of expected deaths, p%) were inversely correlated 
with per capita GDP (r = -0.60), correlated with proportion living in poverty 
(r = 0.66), and modestly correlated with income inequality (r = 0.45). Incidence rate 
ratio for deaths was 1.062 (95% CI, 1.038-1.087) in more versus less vulnerable 
countries. Excess deaths started deviating in the two groups after the first wave. 
Between-country heterogeneity diminished gradually within each group. Less vul- 
nerable countries had mean p% = -0.8% and 0.4% in 0—64 and >65-y-old strata. 
More vulnerable countries had mean p% = 7.0% and 7.2%, respectively. Lower death 
rates were seen in children of age 0-14 y during 2020—2023 versus prepandemic 
years. While the pandemic hit some countries earlier than others, country vulner- 
ability dominated eventually the cumulative impact. Half the analyzed countries 
witnessed no substantial excess deaths versus prepandemic levels, while the others 
suffered major death tolls. 


COVID-19 | mortality | bias | demography 


Estimates of excess deaths are considered to offer an aggregate picture of the impact on 
overall mortality during the COVID-19 crisis (1, 2). Excess deaths capture the composite 
of deaths due to SARS-CoV-2, indirect effects of the pandemic, and the effects of the 
measures taken during the crisis (both positive and negative). Three years after COVID-19 
was proclaimed as a pandemic, one can have a mature picture of the evolution of the 
excess deaths in different countries over time. This is most reliable in countries that have 
adequately complete death registration data and information on the evolution of their 
population age structure, so that proper age-adjustments can be made (3). Performance 
of different countries may have varied markedly. Moreover, excess deaths’ trajectories over 
time have also shown diverse patterns across countries, as different countries witnessed 
more excess deaths during different periods. 

Here, we evaluate excess death trajectories in the period between January 2020 and 
July 2023 in 34 countries with the most reliable data. We evaluate how and why these 
trajectories have diverged across different countries overall and within different age groups. 
We categorize these countries into groups of more and less vulnerable ones. Vulnerability 
is judged based on per capita gross domestic product (GDP), income inequality, and 
proportion of population below the poverty level, since the COVID-19 crisis dispropor- 
tionately affected the poor and the disadvantaged in high-income countries (4-6). 
Examining excess mortality as a function of vulnerability can help gain insights into why 
some countries performed poorly and others did much better during the pandemic. It 
could also offer predictive insights for future pandemic challenges and lessons for pandemic 
preparedness. 


Results 


Cumulative Excess Deaths and Vulnerability Indicators. Table 1 shows the estimated 
cumulative excess deaths (p%, absolute number, and per million) during the entire 
pandemic period and the 3 vulnerability indicators for each of the 34 countries (scatterplot 
in Fig. 1). The Pearson correlation coefficients of p% with GDP per capita, Gini 


PNAS 2023 Vol.120 No.49 e2309557120 


Significance 


Excess deaths during 2020-2023 
reflect direct and indirect effects 
of the COVID-19 pandemic and 
of the measures taken. Data 
from 34 countries with detailed 
death registration and allowing 
adjustment for changes in the 
age structure of the population 
over time show two groups, each 
with very different excess death 
outcomes. The 17 more 
vulnerable countries (those with 
per capita nominal GDP < 
$30,000, Gini > 0.35 for income 
inequality and/or 22.5% of their 
population living in poverty) had 
very high excess deaths 
compared with 2017-2019, while 
the other 17 less vulnerable 
countries had deaths during 
2020-2023 that were comparable 
to 2017-2019. Continuous 
monitoring of excess deaths 
helps understand how country 
vulnerability shapes long-term 
impacts. 
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Table 1. Excess deaths in 2020-2023 (using 2017-2019 average as reference) and vulnerability indicators 


Reported 
Excess COVID-19 
Excess 2021 Excess deaths, deaths per Per capita Gini for 
deaths, population deaths, absolute i million, 2 nominal income Living in 

Country p (%) in millions absolute per million’ August 2023 GDP ($) inequality poverty (%) 
New Zealand =3.56 4.86 -4,619 -950 662 48,781 0.320 No data 
Sweden =3:53 10.16 -11,632 -1,145 2,415 61,029 0.276 0.60 
Denmark =2.98 5.81 -6,245 -1,074 1,509 68,008 0.268 0.40 
Norway -2.05 5.47 -3,161 =578 1,027 89,154 0.263 0.52 
South Korea -1.94 51:31 =22,583 -440 688 34,998 0.331 0.74 
Australia -1.86 25.79 -10,949 -425 872 60,443 0.318 0.74 
Luxembourg -1.80 0.63 =283 -446 1,575 133,590 0.305 0.30 
Iceland =0.39 0.34 =32 =93 542 68,728 0.250 0.04 
Israel 0.48 8.79 829 94 1,433 52,171 0.342 2.47 
Switzerland 0.71 8.72 1,807 207 1,610 91,992 0.316 0.14 
Finland 1.25 5.55 2,547 459 1,818 53,655 0.265 0.10 
Belgium 2.31 11.63 9,224 793 2,955 51,247 0.262 0.20 
Canada 2.66 38.07 25,900 680 1,396 51,988 0.280 0.74 
Germany 2.77 83.90 97,544 1,163 2,086 51,204 0.296 0.50 
France 2.77 65.43 61,191 935 2,568 43,659 0.292 0.11 
Portugal’ 3.08 10.17 12,954 1,274 2,653 24,568 0.310 0.90 
Netherlands 3.35 17.17 19,296 1,124 1,338 57,768 0.304 0.30 
Slovenia‘ 4.00 2.08 3,054 1,469 4,536 29,291 0.246 0.10 
Spain’ 4.15 46.75 64,772 1,386 2,607 30,104 0.320 2.50 
Hungary! 4.57 9.63 21,413 2,223 5,065 18,728 0.286 1.80 
Estonia‘ 4.61 1.33 2,656 2,004 2,188 27,944 0.305 1.30 
UKt 5.04 68.21 114,419 1,678 3,349 46,510 0.355 0.86 
Austria 5.88 9.04 17,605 1,947 2,492 53,638 0.274 0.80 
Italy" 6.11 60.37 140,448 2,327 3,164 35,658 0.330 3.05 
Croatia‘ 6.14 4.08 11,482 2,813 4,480 17,685 0.292 1.80 
Greece’ 6.47 10.37 28,619 2,760 3,593 20,193 0.308 3.40 
Latvia‘ 6.74 1.87 6,793 3,639 3,949 21,148 0.355 2.00 
Slovakia’ 7:27 5.46 14,667 2,686 3,876 21,392 0:222 1.30 
Czechia' 8.30 10.72 33,289 3,104 3,992 26,821 0.248 0.30 
Poland? 8.50 37.80 128,615 3,403 3,165 18,000 0.281 1.20 
Lithuania’ 8.98 2.69 12,314 4,578 3,604 23,723 0.357 1.40 
Chile 10.89 19.21 45,805 2,384 3,207 16,265 0.460 4.40 
USA’ 12.08 332.92 1,220,295 3,665 3,386 70,249 0.375 1.50 
Bulgaria’ 15.76 6.90 59,067 8,565 5,567 12,222 0.408 8.20 


“Excess deaths are calculated from week 1, 2020 until the latest update of available data (Materials and Methods). 
‘Belongs to the group of more vulnerable countries (see Materials and Methods for definitions). 


Downloaded from https://www.pnas.org by 178.204.250.189 on December 27, 2023 from IP address 178.204.250.189. 


2 of 8 


coefficient, and percentage living in poverty were -0.60, 0.45, 
and 0.66. Corresponding Spearman rank correlation coefficients 
were -0.69, 0.29, and 0.62, respectively. The negative binomial 
regression showed that the incidence rate ratio for death was 
1.062 (95% CI 1.038-1.087, P < 0.0001) for more vulnerable 
countries versus less vulnerable countries. Excluding the country 
with the largest population (USA) did not materially change 
the results [1.055 (95% CI, 1.031-1.080, P < 0.0001)]. Each 
of the three economic variables was also significantly associated 
with death risk (P = 0.004 for GDP per capita, P = <0.0001 for 
Gini, P< 0.0001 for percentage living in poverty). The pseudoR” 
for GDP per capita, Gini, and percentage living in poverty was 
0.1653, 0.1712, and 0.1746, respectively. Considering all three 
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variables in a multivariate model did not improve the pseudoR? 
materially (0.1817). 

Across the 34 countries there were 2,097,101 excess deaths, 
58% were accounted by the United States of America alone 
(1,220,295). There were only 176,439 excess deaths in the 17 
countries of the less vulnerable group (500 per million in the 2021 
total population of 352,667,986, people) versus 1,920,662 in the 
17 countries of the more vulnerable group (3,046 per million in 
the total population of 630,541,384 people). 

The United States of America would have had 1.60 million 
fewer deaths if it had the performance of Sweden, 1.07 million 
fewer deaths if it had the performance of Finland, and 0.91 million 
fewer deaths if it had the performance of France. 
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Fig. 1. Scatterplot of excess deaths p% against the per capita GDP also showing countries with substantial income inequality and/or high proportion of 
population living in poverty. The Left panel shows the data for the entire countries’ populations and the other two panels show the strata or nonelderly and 


elderly using a cut-off of 65 y. 


Cumulative Excess Deaths in the Two Vulnerability Groups Over 
Time. Fig. 2 shows country trajectories over time; means and SDs 
for p% for more and less vulnerable countries are in SI Appendix, 
Table S1. Both groups had similar estimates of modest death deficit 
for the first 12 weeks of 2020 (versus the 2017-2019 average) and 
similar SDs. By mid-2020, both groups had reverted from this 
initial death deficit to a picture that resembled 2017-2019 for 
the mean p% (no sizeable excess or deficit), but there was large 
between-country heterogeneity already, especially in the more 
vulnerable group. Six countries had already reached excess deaths 
exceeding 8%. 

After mid-2020, the less vulnerable group of countries maintained 
until July 2023 a very small mean of excess deaths (never exceeding 
1.3%, 0.2% by the latest update). Their between-country heteroge- 
neity also diminished during 2021-2023. The SD at the latest 
update was half the SD of end-2020 in that group. Countries that 
initially had prominent death deficits compared with 2017—2019 
lost much or all of these death deficits by latest follow-up; while 
countries with substantial excess deaths in 2020 or early 2021 
improved subsequently. As of the latest update in 2023, with the 
exception of Austria (p% = 5.9%), all the less vulnerable countries 
had cumulative p% estimates in the narrow range -3.6% to 3.4%. 

Conversely, among more vulnerable countries mean p% reached 
8.2% by end-2020. Between-country heterogeneity diminished 
markedly compared with mid-2020 because although many of 
these countries were spared in the first wave, only Latvia and 
Estonia were still spared by end-2020 (SD of p% fell from 8.2% 
to 4.6%). Mean p% continued to increase during 2021 reaching 
9.9% by mid-2021 and 11.8% by end-2021, while between-country 
heterogeneity remained stable (between 4.4% and 4.8%). During 
2022 and into 2023, mean p% decreased modestly over time and 
some further decrease was seen also in the SD. As of the latest 
update in 2023, with the exception of Portugal (3.1%), all the 
more vulnerable countries had more cumulative excess deaths 
(4.0% to 15.8%) than 16 of the 17 less vulnerable countries. 

There was no correlation between excess deaths p% and per 
capita GDP in 2020; this changed in 2021 and then remained 
steady between -0.55 and -0.62 (SI Appendix, Table S2). Similarly, 
there was no correlation between excess deaths p% and the pro- 
portion of poverty until mid-2020 but this changed in the next 
year and then it remained steady between 0.5 and 0.66. Throughout 
2020-2023 period, there was modest correlation (correlation coef- 
ficients mostly between 0.2 and 0.4) between excess deaths p% 
and income inequality. 


PNAS 2023 Vol.120 No.49 e2309557120 


Sensitivity Analyses Using Week-Specific Trends. In analyses 
using the week-specific trends method, as reported in STMF 
(SI Appendix, Table S3), the excess death estimates for most 
countries tended to be modestly higher than those obtained by 
considering as reference the average of 2017-2019 (SI Appendix, 
Fig. S1). Nevertheless, the overall patterns (difference between less 
and more vulnerable, correlations with economic indicators, and 
evolution over time) remained largely unchanged. Countries in 
the two groups were separated by a cut-off of p% = 9.1%, except 
for Netherlands (13.7%) and Austria (9.4%) having higher values 
in the less vulnerable group and Spain (8.1%), Lithuania (7.5%), 
and Italy (7%) having lower values in the more vulnerable group. 


Age Groups. Across all 34 countries, 28% of total excess deaths 
were in the nonelderly (0—64 y) and 72% in the elderly (65 and 
over) strata (Table 2). For the nonelderly, at the latest update, 
excess deaths means were -0.8% and 7.0% in the less and more 
vulnerable groups, respectively. For the elderly, respective means 
were 0.4% and 7.2%. In the 17 less vulnerable countries, there 
were 13,472 and 162,967 excess deaths in the nonelderly and 
nonelderly strata, respectively; nonelderly excess deaths accounted 
for 7.6% of total excess deaths. In the 17 more vulnerable 
countries, respective numbers were 577,281 and 1,343,381; 
nonelderly excess deaths accounted for 30% of total excess deaths. 

Most countries had slightly lower p% in the nonelderly and 
elderly strata (by 1.8%), with some notable exceptions (Table 2 
and Fig. 1). Canada, United Kingdom, United States of America, 
and Chile had much higher p% in the nonelderly than the elderly 
(difference over 7.2%). The United States of America had the 
highest excess deaths than any country among the nonelderly 
(p% = 17.9%), while Canada had high excess deaths in the 
nonelderly (p% = 9.5%), but no excess deaths in the elderly 
(p% = 1.1%). Higher estimates in the nonelderly rather than 
elderly were manifest in these countries from the first pandemic 
months and persisted subsequently. Conversely, Slovenia, France, 
and Poland performed much better in the nonelderly than in the 
elderly (difference over 6.9%); France and Slovenia had death 
deficits in nonelderly. 

For children (0-14 y), most countries had lower death rates 
during 2020-2023 versus 2017-2019. Numerically, the largest 
death deficits came actually from more vulnerable countries [USA 
(-3,456), Chile (-1,318), Italy (-1,094), and Poland (-996)]. 

SI Appendix, Table S4 provides estimates of p%, absolute 


excess deaths and per million excess deaths in each country for 
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Estimated cumulative excess deaths at different timepoints during 2020-2023 in 34 countries grouped into more vulnerable (upper group) and less 


vulnerable (lower group). The color scale corresponds to excess deaths increasing from green to red. 


each of the 2 nonelderly age strata (I-14, 15-42) and SI Appendix, 
Table S5 does the same for each of the 3 elderly strata (65-74, 
75-84, and 85+). SI Appendix, Table S6 provides standardized 
excess death rates based on the 2013 revision of the European 
Standard Population. 


Discussion 


Analysis of data during 2020-2023 from 34 countries with reliable 
mortality data and information on population age structure shows 
that half of these countries had minimal excess deaths or even death 
deficits during 2020-2023 versus the 3 pre-pandemic years, while 
the other half had substantial excess deaths. Performance of different 
countries could be grouped based on vulnerability indicators per- 
taining to their wealth, income inequality, and poverty. Excess 
deaths during 2020-2023 were strongly inversely correlated with 
per capita GDB strongly correlated with proportion of population 
living in poverty and modestly correlated with income inequality. 
‘These correlational patterns were not seen during the first wave, but 
they became manifest and persistent subsequently. 
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Some less vulnerable countries transiently had notable excess deaths 
during 2020 and/or 2021. However, by mid-2023, all except Austria 
(and Netherlands in a week-specific trend sensitivity analysis, 
SI Appendix, Fig. S2) had no major excess deaths considering the 
entire 2020-2023 period. Perhaps transient peaks of excess deaths in 
these countries were to a large extent due to the demise of frail elderly 
individuals, e.g., in long-term care facilities, with limited life expec- 
tancy (7). For individuals with life expectancy less than 2 to 3 y, 
premature deaths in the early pandemic, contribute no excess deaths 
when the whole 2020-2023 period is considered (7, 8). The group 
of less vulnerable countries have high income, no prominent income 
inequalities and no major proportion of their population living in 
poverty. All of them also have practically universal health coverage. 
Thus, they may have fewer excess deaths among disadvantaged people. 
‘These countries also had more means to mobilize sufficient healthcare 
and public health resources. The exact ingredients of their success 
cannot be deciphered from ecological, country-level data. Regardless, 
despite having many elderly, these countries managed to pass the crisis 
with no more deaths than during the recent prepandemic years, after 
adjusting for age structure. 
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Table 2. Excess deaths from January 2020 and up to latest update (week 29, 2023) per age group’ 


Country 0-64 y, absolute 0-64y,p% 65y and over, absolute 65yandover,p% 0-14y, absolute 0-14y, p% 
Sweden -2,388 -7.02 -9,244 -3.1 -80 -5.8 
Denmark -2,022 -6.95 -4,223 -2.3 -48 -4.9 
France -14,323 -4.2 75,514 4.0 -618 -4.9 
Luxembourg -99 -3.7 -184 -1.4 NR NR 
Finland -944 -3.4 3,491 2.0 -45 -7.8 
Slovenia -360 -3.1 3,414 53 -18 -8.3 
Norway =621 =3,1 -2,540 =1.9 -69 =10:3 
Israel -809 =2.7 1,638 1A -417 -14.4 
New Zealand -548 -2.3 -4,071 -3.9 -41 -3.4 
Belgium -1,221 -2.1 10,445 3.1 -460 -21.8 
Australia -2,004 -2.0 -8,945 -1.8 -562 =12.1 
Switzerland -421 -1.4 2,228 1.0 -48 -3.4 
South Korea 1,338 0.6 -23,921 -2.6 -938 -19.2 
Hungary 668 0.7 20,745 5.6 -150 -8.8 
Netherlands 1,251 1:7 18,045 3.6 4 0.2 
Croatia 739 25 10,743 6.8 -41 -5.7 
Poland 10,433 3.1 118,182 10.0 -996 -13.2 
Germany 16,621 3.5 80,923 2.7 -750 -5.8 
Iceland 47 3.5 -79 -1.1 NR NR 
Portugal 2,443 4.3 10,511 2.9 -107 -8.0 
Czechia 3,028 4.6 30,261 9.0 =272 -17.0 
Estonia 572 5.4 2,084 4.4 =12 -6.5 
Slovakia 2,546 5.5 12,121 7.8 -60 =4;1 
Italy 13,457 5.7 126,991 6.2 -1,094 -17.2 
Spain 11,819 5.7 52,953 3.9 -221 -4.2 
Austria 2,568 6.2 15,037 5.8 -118 -9.3 
Latvia 1,539 6.9 5,254 6.7 -97 -24.8 
Greece 4,536 8.3 24,083 6.2 -170 =i 
Canada’ 17,047 9.5 8,853 1.1 1,072 14.5 
Lithuania 2,924 9.7 9,390 8.8 -108 -21.2 
UK 39,740 11.8 74,679 3.9 -823 -5.9 
Bulgaria 10,795 14.3 48,272 16.1 -318 -16.5 
Chile 17,038 16.3 28,767 9.1 -1,318 -18.7 
USA 455,364 17.9 764,931 10.1 -3,456 -3.3 


*Countries are listed in increasing p% among the 0-64-y-old age stratum. 
‘For Canada, the 0-14-y-old age band given in stmf is extrapolated from the 0-44-y age band in the original data https://www150.statcan.gc.ca/n1/tbl/csv/13100768-eng.zip. See https:// 
www.mortality.org/File/GetDocument/Public/STMF/DOC/STMFNote.pdf for stmf methods. 
NR: not reliable (given the very small population and extremely small number of deaths of children, excess deaths can be markedly influenced by single anomalies, e.g., in Iceland there 
were only 6 reported deaths in that age stratum in 2019 (less than a third versus 2017) and in Luxembourg there were 76 in 2022 (more than double versus the 2017-2019 average). 


Conversely, more vulnerable countries fared poorly during 


2020-2023. Their excess death estimates have decreased modestly 
after peaking by end-2021. This may reflect a relative contribution 
of the same phenomenon of some early deaths among frail elderly 
people with limited life expectancy, especially in Slovenia, Slovakia, 
Poland, Spain, and United Kingdom, where decreases in cumu- 
lative excess deaths were prominent over time. It may also reflect 
the more limited fatality of Omicron waves in 2022-2023 (9), 
especially given almost ubiquitous prior infections and high rates 
of vaccination (10). However, cumulative excess death estimates 
remain high (and occasionally very high) in the more vulnerable 
countries. Some countries already reached high excess levels in 
early 2020, while others succumbed during the second wave. 
Again, the exact ingredients of their failure cannot be pinpointed 
in detail. However, the common denominator for these countries 
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was their less robust economies and/or large share of poor, disad- 
vantaged people. COVID-19 was a crisis of inequalities and many 
measures taken may even have fostered worsening inequalities. 
The United States of America is a striking case, with extremely 
high cumulative excess death rates despite high per capita GDP. 
USA income inequality is high, many people live in poverty, and 
many lack health insurance coverage (among working-age people 
15% in 2019 and 12% in 2022, among the entire population 
10% in 2019 and 8% in 2022) (11). Up to 27-fold differences in 
COVID-19 death rates in strata defined by race, gender, and 
educational achievement have been described in USA (6). Area 
deprivation may be a risk factor for COVID-19 mortality, espe- 
cially among minorities (12, 13). Moreover, many excess deaths 
apparently are due to sharply increasing deaths due to overdose 
(14) and deaths due to suboptimal healthcare access during 
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2020-2023 (15). Such noninfectious causes may largely explain 
the exceptionally high excess death p% estimate among USA 
nonelderly people. 

Income inequalities have been strongly associated with COVID-19 
deaths also in other vulnerable countries, e.g., in Spain (16, 17) and 
Chile (18). The pandemic widened all-cause preexisting mortality 
gaps according to income even in some less vulnerable countries, e.g., 
Netherlands (19) and Sweden (20). In a vicious circle, the pandemic 
caused more poverty worldwide (21). 

‘The gap between the more and less vulnerable group of countries 
was overall similar among the elderly and the nonelderly. However, 
Canada, United Kingdom, United States of America, and Chile 
had very poor performance specifically among the nonelderly, while 
Slovenia, France, and Poland showed the opposite pattern. 
Population obesity rates are markedly different in these two groups 
of countries (28-36% versus 20-23%) and additional differences 
in background health status of nonelderly individuals may explain 
these patterns, e.g., overdose and poor healthcare access in the 
United States of America—as discussed above. The United States 
of America, Chile, and the United Kingdom were categorized here 
as more vulnerable countries, while Canada was categorized as less 
vulnerable. However, the classification is not absolute. Canada also 
has a modest proportion of people living in poverty (0.74%), and 
data from Ontario suggest that low income and low educational 
attainment increased risks of COVID-19 deaths (22). In British 
Columbia, death rates from overdose more than doubled during 
the pandemic versus 2019 (23). 

For children, most countries had fewer deaths during 2020- 
2023 compared with the 3 prepandemic years. This may be due 
to the exceptionally low infection fatality rate of SARS-CoV-2 
among children (24), plus the almost complete disappearance of 
influenza for 2 y. More granular data should be examined to see 
whether death deficits extend also to adolescents and very young 
adults. 

Our study has some limitations. First, excess death estimates 
depend on modeling assumptions that may affect the absolute 
values (25, 26). In a sensitivity analysis considering the SMTF 
week-specific trends with decreasing mortality rates over the years, 
the estimated excess deaths were higher than what we calculated 
using the 2017-2019 average as reference. However, even if mor- 
tality rates had decreased in many countries over time in the past, 
there is no guarantee that decreases should continue, especially 
among the increasingly elderly and frail populations of high-income 
countries (27). In essence, the analysis including linear trends 
based on past years assumes that a continued 2—5% decrease in 
death rate should continue to materialize every year in most of 
these high-income countries. When there is no increase in deaths, 
this is then counted as substantial excess mortality when calibrated 
against that optimistic extrapolation of continuous constant 
improvement. Furthermore, extrapolating the trend 4 y out from 
the end of the reference period may also be problematic for any 
trend analysis. Moreover, we have previously shown that different 
modeling assumptions do not affect the comparative performance 
of different countries (27): Poor performers are consistently poor 
performers regardless of the exact modeling and this was largely 
borne out also in the current time-trend sensitivity analysis. 
Therefore, the group-level contrast between less and more vulner- 
able countries is probably robust to modeling choices. Of note, 
one should avoid making inferences about the relative ranking of 
countries in pairwise comparisons. The p% for each country are 
estimates and they carry large uncertainty, not only because of 
statistical chance error (which may be modeled by a Poisson or 
negative binomial process for expected deaths) but also additional 
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uncertainty due to potential errors in the underlying data and 
multiplicity of options in the modeling assumptions. Instead, one 
should focus on the big picture of the relative performance of the 
two large groups of less and more vulnerable countries. 

Second, some data may have imperfections, e.g., some missing 
deaths in the last few weeks of the covered period, or inaccuracies 
in inferred data for the population structure in each age bin and 
country over time. However, the data for the 34 chosen countries 
are likely to be the most reliable. Extrapolations to countries with 
less reliable death and age structure data are precarious (28). 
Socioeconomic factors may have had a major impact on mortality 
during 2020-2023 also in other countries, perhaps even to a larger 
extent than in the 34 countries analyzed here (29-31). Their 
impact might have been even larger for non-COVID-19 deaths 
(29). Conversely, most countries not analyzed here have much 
younger populations and few frail elderly in long-term care facil- 
ities than the 34 analyzed countries; thus, deaths directly due to 
SARS-CoV-2 are expected to have been fewer (24, 32). 

Third, we used three high-level economic indicators and these 
should be seen as surrogate of a mixture of many socioeconomic and 
other factors that operate at individual, household, work, community, 
and societal levels. Poor outcomes during the COVID-19 crisis are 
probably determined by a large web of interrelated factors. 

Acknowledging these caveats, the country-level analyses offer 
a picture of how countries succeeded or failed in dealing with this 
major crisis. Half the countries went through the 3 crisis years 
without witnessing substantial excess deaths versus their 2017— 
2019 prepandemic levels, while others suffered major death tolls. 
Continuous monitoring of excess death patterns may be useful, 
given that the impact of the COVID-19 crisis and several measures 
taken may continue for several years. 


Materials and Methods 


Excess Death Calculations, Eligible Countries, and Mortality Data. We 
performed excess death calculations using methods similar to those we have 
used in our previously published work for 2020-2021 (3). Briefly, we used the 
average of years 2017-2019 as baseline and included age-adjustment in the 
calculations, considering data in 5 age strata (0-14, 15-64, 65-75, 75-85, and 
>85 y old). For each age stratum, we obtained the average mortality, the number 
of deaths per million for the population of the specific age stratum and estimated 
the expected deaths during each week of the pandemic periods of interest, cor- 
recting for population size in the specific age stratum and summing expected 
deaths across population strata. 

We used the Human Mortality Database (https://www.mortality.org) (33), 
specifically the Short Term Mortality Fluctuations (STMF) file (https://www.mortality. 
org/File/GetDocument/Public/STMF/Outputs/stmf.csv) that includes weekly data. 
Only countries that have excellent death registration and include data with weekly 
deaths in the STMF during the period from 2017 until at least the second half of 
2022 were considered, to avoid spurious changes in recorded deaths over time due 
to changes in death registration. STMF had data from 34 such eligible countries 
when we downloaded the database (August 2023). Data covered up to at least 
week 20 (at most up to week 29, 2023) for 32 countries, up to week 17 for Australia, 
and up to week 8 for Canada. 

We calculated weekly excess deaths and focused on cumulative excess deaths 
up to the end of each week starting from the first week of 2020. Calculations 
proceeded until the latest week with available death data in each country. The 
Human Mortality database uses ISO weeks that end on Sunday, e.g., the first 
week of 2020 ends on Sunday, 5 January 2020. We expressed the excess death 
impact as percentage above expected deaths, p%. 

All analyses use 2017-2019 as the reference period to calculate expected deaths. 
We have previously shown that relative country performance tends to be similar with 
shorter and longer windows for the reference period (27), even if exact absolute 
excess deaths may vary with different modeling options (27). There is also debate 
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in the literature about if and how linear (or other) trends should be considered in 
the modeling, with some researchers arguing in favor of considering such trends 
(34, 35).We have previously shown that once changes in the age pyramid are taken 
into account, extra consideration of trends may be unnecessary or even unwanted 
(e.g., if it extrapolates that future mortality trend patterns should continue in the 
future, even if they are saturated (3, 27). Nevertheless, for comparison, we also 
provide excess death estimates for the whole 2020-2023 period using the week- 
specific trends method that is already implemented in STMF (7, 36); all available 
consecutive years in the window 2000-2019 are considered for the reference period 
and week-specific linear trends are considered in this option. 


Classification of Countries in Vulnerability Groups. We aimed to classify the 
34 countries into two groups of 17 countries each using a parsimonious set of 
economic variables to separate a more vulnerable group from a less vulnerable 
group. We reasoned that per capita nominal GDP should be a most straight- 
forward classifier to consider. However, we also reasoned that one should also 
allow for high vulnerability when, despite a high per capita nominal GPD, there 
is large income inequality and/or there are many people who are very poor and 
disadvantaged and thus would not be shielded by the services that can be offered 
otherwise to the majority of citizens. Indeed, it is well known that the Preston curve 
(the relationship between income and life expectancy) (37) is almost flat after a 
certain income (i.e., further increases in income do not improve life expectancy) 
and that disadvantaged segments of the population can have much higher mor- 
tality, even in wealthy countries. Therefore, we defined a priori the more vulnera- 
ble countries’ group to include countries with per capita nominal GDP < $30,000, 
large income inequality (Gini > 0.35), and/or sizeable percentage (=2.5%) of 
population in poverty [defined as living with <$5.50 per day(2011 PPP values)]. 
The group of less vulnerable countries had none of these vulnerability features. 
These three indicators may correlate with many other indicators of vulnerability 
that span the individual, household, work, community, and societal levels. They 
may also offer surrogacy for lack of health care resources for poor, disadvantaged 
segments of the population. 

Information on these three indicators was obtained from World Bank for per cap- 
ita nominal GDP [2021 estimates (38); OECD for Gini of income inequality (2020 or 
more recent estimates, “Income Distribution Database.” OECD.org). Measure: Gini 
(disposable income, post taxes and transfers), coupled with data from the World 
Bank for countries that are not OECD members (39, 40); and Wikipedia (2019 for 
most recent data of proportion of population living in poverty) (41). 


1. M.V. Kiang, R.A. Irizarry, C. O. Buckee, S. Balsari, Every body counts: Measuring mortality from the 
COVID-19 pandemic. Ann. Intern. Med. 173, 1004-1007 (2020). 

2. W. Msemburi et al., The WHO estimates of excess mortality associated with the COVID-19 pandemic. 
Nature 613, 130-137 (2023). 

3. M. Levitt, F. Zonta, J. P. A. loannidis, Comparison of pandemic excess mortality in 2020-2021 across 
different empirical calculations. Environ. Res. 213, 113754 (2022). 

4. M.S. Shiels et al., Racial and ethnic disparities in excess deaths during the COVID-19 pandemic, 
March to December 2020. Ann. Intern. Med. 174, 1693-1699 (2021). 

5. Y. Commodore-Mensah, L. A. Cooper, Reversing the tide of racial and ethnic disparities in excess 
deaths during the COVID-19 pandemic. Ann. Intern. Med. 174, 1755-1756 (2021). 

6. E.B. Pathak, J. M. Menard, R. B. Garcia, J. L. Salemi, Joint effects of socioeconomic position, 
race/ethnicity, and gender on COVID-19 mortality among working-age adults in the United States. 
Int. J. Environ. Res. Public Health 19, 5479 (2022). 

7. M. Ballin etal., Time-varying risk of death after SARS-CoV-2 infection in Swedish long-term care 
facility residents: A matched cohort study. BMJ Open. 12, e066258 (2022). 

8. F. Espuny Pujol, R. Hancock, M. Morciano, Trends in survival of older care home residents in 
England: A 10-year multi-cohort study. Soc. Sci. Med. 282, 113883 (2021). 

9. N.U. Friis et al., COVID-19 mortality attenuated during widespread Omicron transmission, Denmark, 
2020 to 2022. Eur. Surveill. 28, 2200547 (2023). 

10. S. Pilz, J. P.A. Ioannidis, Does natural and hybrid immunity obviate the need for frequent vaccine 
boosters against SARS-CoV-2 in the endemic phase? Eur. J. Clin. Invest. 53, 13906 (2023). 

1. R.A. Cohen, A. E. Cha, Health insurance coverage: Early release of estimates from the National 
Health Interview Survey (2022). https://www.cdc.gov/nchs/data/nhis/earlyrelease/insur202305_1. 
pdf. Accessed 8 August 2023. 

2. S.Magesh et al., Disparities in COVID-19 outcomes by race, ethnicity, and socioeconomic status: 
Asystematic-review and meta-analysis. JAMA Network Open 4, 2134147 (2021). 

3. E.C. Matthay et al., Occupation and educational attainment characteristics associated with 
COVID-19 mortality by race and ethnicity in California. JAMA Network Open 5, e228406 (2022). 

14. https://www.cdc.gov/nchs/pressroom/nchs_press_releases/2022/202205.htm. Accessed 26 March 2023. 

5. G.A. Roth, M. Vaduganathan, G. A. Mensah, Impact of the COVID-19 pandemic on cardiovascular 
health in 2020: JACC state-of-the-art review. J. Am. Coll. Cardiol. 80, 631-640 (2022). 

16. J. Politi et al., Epidemiological characteristics and outcomes of COVID-19 cases: Mortality 

inequalities by socio-economic status, Barcelona, Spain, 24 February to 4 May 2020. Eur. Surveill. 
26, 2001138 (2021). 


PNAS 2023 Vol.120 No.49 e2309557120 


The magnitude of the differences across more and less vulnerable countries 
can be best appreciated if one estimates how many fewer deaths a more vulner- 
able country would have had if it had the same performance as a less vulnerable 
country. To obtain this number, the absolute excess deaths of the more vulnerable 
country are multiplied by the difference in p% between the compared countries. 
IIlustratively, we present such absolute numbers of deaths that could have been 
saved for USA, if it had the performance of Sweden, Finland, or France. 


Analyses. Cumulative excess deaths for each country were analyzed for each 
week starting with week 1, 2020 and until the most recent available data. For 
the few countries with missing data on the last few weeks in this interval, the last 
available value was carried forward. Main milestones were: week 12 (22 March 
2020, at the start of the impact of the pandemic wave in Europe and USA), mid- 
2020, end-2020, mid-2021, end-2021, mid-2022, end-2022, mid-2023, and 
latest available (up to at most week 29, 2023). 

We examined Pearson and Spearman correlation coefficients for the correla- 
tion between p% at each main milestone and the three economic variables (GDP 
per capita, Gini coefficient, percentage of population living in poverty) across 
all 34 countries. We also performed a negative binomial regression (to allow 
for overdispersion) to examine the impact of vulnerable status and each of the 
three economic indicators. The negative binomial regression used population 
denominations as exposure and also considered the expected rates of deaths. 

We performed similar analyses splitting excess death estimates in nonelderly 
(0-64 y, i.e. combining the 0-14 and 15-64 y old strata) and elderly people (65 y 
old and above, i.e., combining the 65-74, 75-84, and over 85 y old strata). We 
also examined children in particular (0-14 y). Estimates of p%, absolute and per 
million excess deaths for all granular age strata are provided in supplements so 
as to allow any type of standardized comparisons (42) by interested researchers. 


Data, Materials, and Software Availability. All study data are included in the 
article and/or S/ Appendix. 
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